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This Colloquium reviews statistical models for money, wealth, and income distributions devel- 
oped in the econophysics literature since the late 1990s. By analogy with the Boltzmann-Gibbs 
distribution of energy in physics, it is shown that the probability distribution of money is 
exponential for certain classes of models with interacting economic agents. Alternative scenarios 
are also reviewed. Data analysis of the empirical distributions of wealth and income reveals a 
two-class distribution. The majority of the population belongs to the lower class, characterized 
by the exponential ("thermal") distribution, whereas a small fraction of the population in the 
upper class is characterized by the power-law ("superthermal") distribution. The lower part is 
very stable, stationary in time, whereas the upper part is highly dynamical and out of equilibrium. 

"Money, it's a gas." Pink Floyd, Dark Side of the Moon 
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I. HISTORICAL INTRODUCTION 

This Colloquium article is based on the lectures that 
one of us (V.M.Y.) has frequently given during the 
last nine years, when econophysics became a popular 
subject. Econophysics is a new interdisciplinary re- 
search field applying methods of statistical physics to 
problems in economics and finance. The term "econo- 
physics" was first introduced by the theoretical physi- 
cist Eugene Stanley in 1995 at the conference Dynam- 
ics of Complex Systems, which was held in Kolkata as 
a satell ite meeting to t he STATPHYS- 19 confer ence in 
China (|Carbone et all 120071: IChakrabartil I2005D . The 



term appeared first bv lStanlev et all ( 19961 ) in the pro- 
ceedings of the Kolkata conference. The paper pre- 
sented a manifesto of the new field, arguing that "be- 
havior of large numbers of humans (as measured, e.g., 
by economic indices) might conform to analogs of the 
scaling laws that have proved useful in describing sys- 
tems composed of la rge numbers of inanimate objects" 
( Stanlev et a/.l . [19961 ). Soon the first econophysics confer- 
ences were organized: International Workshop on Econo- 
physics, Budapest, 1997 and International Workshop 
on Econophysics and Statistical Finance, Palermo, 1998 
( Carbone et all. l2007l). and the book An Introduction to 
Econophysics by Mantegna and StanlevI ( 1999( ) was pub- 
lished. 

The term econophysics was introduced by analogy with 
similar terms, such as astrophysics, geophysics, and bio- 
physics, which describe applications of physics to differ- 
ent fields. Particularly important is the parallel with 
biophysics, which studies living organisms, but they still 
obey the laws of physics. Econophysics does not literally 
apply the laws of physics, such as Newton's laws or quan- 
tum mechanics, to humans. It uses mathematical meth- 
ods developed in statistical physics to study statistical 
properties of complex economic systems consisting of a 
large number of humans. As such, it may be considered 
as a branch of applied theory of probabilities. However, 
statistical physics is distinctly different from mathemat- 
ical statistics in its focus, methods, and results. 

Originating from physics as a quantitative science, 
econophysics emphasizes quantitative analysis of large 
amounts of economic and financial data, which became 
increasingly available with the introduction of comput- 
ers and the Internet. Econophysics distances itself from 
the verbose, narrative, and ideological style of political 
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economy and is closer to econometrics in its focus. Study- 
ing mathematical models of a large number of interacting 
economic agents, econophysics has much common ground 
with the agent-based modeling and simulation. Corre- 
spondingly, it distances itself from the representative- 
agent approach of traditional economics, which, by def- 
inition, ignores statistical and heterogeneous aspects of 
the economy. 



Another direction related to econophysics has been ad- 
vocated by the theoretical physicist Serg e Galam sinc e 
early 1980 under the name of sociophysics (|GalamV2004l). 
with the first appearance of the term by [Galam et al. 
(|1982| ). It echoes the term "physique sociale" pro- 
posed in the nineteenth century by Auguste Comte, the 
founder of sociology. Unlike econophysics, the term "so- 
ciophysics" did not catch on when first introduced, but 
it is coming back with the popularity of e conophysics 
and a.ctive support from sorne phy sicists ( Schweitzer] . 
120031: IStaufferl . l2004l : IWeidlichl . |2000D . While the prin- 
ciples of both fields have much in common, econophysics 
focuses on the narrower subject of economic behavior 
of humans, where more quantitative data is available, 
whereas sociophysics studies a broader range of social 
issues. The boundary between econophysics and socio- 
physi cs is not sharp, and the two fields enjoy a good ra p- 
port ( Chakrabarti. Chakraborti. and Chatteried , |2006( ). 



Historically, statistical mechanics was developed in the 
second half of the nineteenth century by James Clerk 
Maxwell, Ludwig Boltzmann, and Josiah Willard Gibbs. 
These physicists believed in the existence of atoms and 
developed mathematical methods for describing their sta- 
tistical properties. There are interesting connections be- 
tween the development of statistical physics and statistics 
of social phenomena, which were recently highlighted by 
the science journalist Philip BaU (2002, 2004). 

Collection and study of "social numbers" , such as the 
rates of death, birth, and marriage, has been 
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progressively since the seventeenth century (|Ball[ 
Ch. 3). The term "statistics" was introduced in the eigh- 
teenth century to denote these studies dealing with the 
civil "states" , and its practitioners were called "statists" . 
Popularization of social statistics in the nineteenth cen- 
tury is particularly accredited to the Belgian astronomer 
Adolphe Quetelet. Before the 1850s, statistics was con- 
sidered an empirical arm of political economy, but then 
it started to transform into a general method of quanti- 
tative analysis suitable for all disciplines. It stimulated 
physicists to develop statistical mechanics in the second 
half of the nineteenth century. 

Rudolf Clausius started development of the kinetic the- 
ory of gases, but it was James Clerk Maxwell who made 
a decisive step of deriving the probability distribution of 
velocities of molecules in a gas. Historical studies show 
(|Ball, 200-1, Ch. 3) that, in developing statistical mechan- 
ics, Maxwell was strongly influenced and encouraged by 
the widespread popularity of social statistics at the time 



('Gillispiel. [T963h .^ This approach was further developed 
by Lud wig Boltzm ann, who was very explicit about its 
origins (|Bal]| . |2004 p. 69): 



"The molecules are like individuals, . . . and 
the properties of gases only remain unaltered, 
because the number of these molecules, which 
on the average have a given state, is con- 
stant." 



In his book Populare Schrifen. iBoltzmannI ( 1905| ) praises 
Josiah Willard Gibbs for systematic development of sta- 
tistical mechanics. Then, Boltzmann says:^ 

"This opens a broad perspective, if we do not 
only think of mechanical objects. Let's con- 
sider to apply this method to the statistics of 
living beings, society, sociology and so forth." 

It is worth noting that many now-famous economists 
were originally educated in physics and engineering. Vil- 
fredo Pareto earned a degree in mathematical sciences 
and a doctorate in engineering. Working as a civil engi- 
neer, he collected statistics demonstrating that distribu- 
tions of incom e and wealth in a society follow a power law 
(lParetolfl897l) . He later became a professor of economics 
at Lausanne, where he replaced Leon Walras, also an en- 
gineer by education. The influential American economist 
Irving Fisher was a student of Gibbs. However, most 
of the mathematical apparatus transferred to economics 
from physics was that of New tonian mechanics and clas- 
sical thermodynamics ( Mirow ski. 1989; Smith and Folc;^, 
l2008l) . It culminated in the neoclassical concept of mech- 
anistic equilibrium where the "forces" of supply and de- 
mand balance each other. The more general concept 
of statistical equilibrium largely eluded mainstream eco- 
nomics. 

With time, both physics and economics became more 
formal and rigid in their specializations, and the social 
origin of statistical physics was for gotten. Th e situation 
is weh summarized by Phihp BaU jBalll . |2004 p. 69): 



"Today physicists regard the application of 
statistical mechanics to social phenomena as 
a new and risky venture. Few, it seems, re- 
call how the process originated the other way 
around, in the days when physical science 
and social science were the twin siblings of a 
mechanistic philosophy and when it was not 
in the least disreputable to invoke the habits 
of people to explain the habits of inanimate 
particles." 

Some physicists and economists attempted to connect 
the t wo disciplines during the twentieth century. Fred- 
erick [Sodd^ (I1933D . the Nobel Prize winner in chemistry 



^ V.M.Y. is g rateful to Stephen G. Brush for this reference. 
Cited from IBoltzmannI l|2006h . V.M.Y. is grateful to Michael 
E. Fisher for this quote. 
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for his work on radioactivity, published the book Wealth, 
Virtual Wealth and Debt, where he argued that the real 
wealth is derived from the energy use in transforming raw 
materials into goods and services, and not from monetary 
transactions. He also warned about dangers of exces- 
sive debt and related "virtual wealth" , thus anticipating 
the Great Depression. His ideas were largely ignored at 
the time, but resonate todav (,Defilla,, 2007) . The theo- 
retical physicist Ettore iMaioranal |l9421) argued in favor 
of applying the laws of statistical physics to social phe- 
nomena in a paper published after his mysterious disap- 
pearance. The statistical physicist Elliott MontroU co- 
authored the book In troduction to Qua ntitative A spects 
of Social Pheno rnena (iMontroll a nd Ba dgeij . 1974 ) Sev- 
eral economist s (|Blumel . ll993HDurlauill997[lFolevl . [l99l 
iFollmeij . I1974J ) applied statistic a.1 physics to economic 
proble ms. The mathematicians iFarioun and Machoved 
(I1983D argued that many paradoxes in classical polit- 
ical economy can be resolved if one adopts a prob- 
abilistic approach. An early attempt to bring to- 
gether the leading theoretical physicists and economists 
at the Santa Fe Institute was not entirely successful 
( Anderson. Arrow, and Pinei . Il988f ). However, by the 
late 1990s, the attempts to apply statistical physics to 
social phenomena finally coalesced into the robust move- 
ments of econophysics and sociophysics. 

Current standing of econophysics within the physics 
and economics communities is mixed. Although an en- 
try on econophysics has a ppeared in the New Palgrave 
Dictionary of Economics ( Rosseii l2008afl , it is fair to 
say that econophysics has not been accepted yet by 
mainstream economics. Nevertheless, a number of open- 
minded, nontraditional economists have joined this move- 
ment, and the number is growing. Under these cir- 
cumstances, econophysicists have most of their papers 
published in physics journals. The journal Physica A: 
Statistical Mechanics and its Applications has emerged 
as the leader in econophysics publications and has even 
attracted submissions from some bona fide economists. 
Gradually, reputable economics jou rnals are also start 
ing to pub lish econophysics papers (Gabaix et at, 200| 
ILux and Sornette . 2002; Wyart and Bouchaud, 200 
The mainstream physics community is generally sym- 
pathetic to econophysics, although it is not uncom- 
mon for econophysics papers to be rejected by Phys- 
ical Review Letters on the grounds that "it is not 
physics" . There is a PACS number for econophysics, 
and Physical Review E has published many papers on 
this subject. There are regular conferences on econo- 
physics, such as Applications of Physics in Financial 
Analysis (sponsored by the European Physical Soci- 
ety) , Nikkei Econophysics Symposium, Econophysic s Col- 
loquium, and Econophys-Kolkata CChak rabartij, l2005t 
IChatteriee. Yarlagadda. and Chakrabartii[2005[ ). Econo- 
physics sessions are included in the annual meetings of 
physical societies and statistical physics conferences. The 
overlap with economists is the strongest in the field of 
agent-based simulation. Not surprisingly, the conference 



series WEHIA/ESHIA, which deals with heterogeneous 
interacting agents, regularly includes sessions on econo- 
ph ysics. More information can be fo u nd in the reviews 
bv jFarmer. Shubik. and SmithI (|2005l) : ISam anido u et oil 
(|2007l ) and on the Web portal Econophysics Forum 



http : //www . unif r . ch/ econophysics/ 



II. STATISTICAL MECHANICS OF MONEY 
DISTRIBUTION 

When modern econophysics started in the middle of 
1990s, its attention was primarily focused on analysis of 
financial markets. Soon after, another direction, closer 
to economics than finance, has emerged. It studies the 
probability distributions of money, wealth, and income 
in a society and overlaps with the long-standing line of 
research in economics studying inequal ity in a society.^ 
Many papers in the economi c litera ture ( Champernownd . 
I1953I : iGibratl . Il93ll: iKaleckil Il945^ use a stochastic pro- 
cess to describe dynamics of individual wealth or income 
and to derive their probability distributions. One might 
call this a one-body approach, because wealth and 
income fluctuations are considered independently for 
each economic agent. Inspired by Boltzmann's kinetic 
theory of collisions in gases, econophysicists introduced 
an alternative, two-body approach, where agents per- 
form pairwise economic transactions and transfer money 
from one agent to another. Actually, this approach 
was pioneered by the sociologist John lAngle (Ji)8^ 
11991 119931 11996L I2002D already in the 1980s. However, 
his work was largely unknown until it was brought 
to the attention of econophysicists by the economist 
Thomas ILuxI (j2005l ) . Now, A ngle's work is widely 



cited in econophvsics literature (lAngk . 20061). Mean- 
while, the physicists llspolatov. Krapivskv. and Rednej 
( 19981 ) independently introduced a statistical model 



of pairwise money transfer between economic agents, 
which is equivalent to the mo del of Angle. Soon, 
three influential papers by Boucha ud and MezardI 



(I2OOOI): IChakraborti and Chak rabartl (|2000l ); 
iDragulescu and Yakovenkol ( 2000t) appeared and 



generated an expanding wave of follow-up publications. 
For pedagogical reasons, we start reviewing this subject 
with the simplest ver sion of the pairwise money t ransfe r 
models presented in IDragulescu and Yakovenkol ( 2000t ). 
This model is the most closely related to the traditional 
statistical mechanics, which we briefly review first. Then 
we discuss the other models mentioned above, as well as 
numerous follow-up papers. 

Interestingly, the study of pairwise money trans- 
fer and the resulting statistical distribution of money 



^ See, e.g., [Atkinson and BourEuignonl 

Atkinson and Pikettvl ll2007j): Ichampernownj 



iJOOO") 
i l953il 



Champernowne and Cowelll ill 9981'): IGibratl lll93lf): lKakwa,nil 
| |1980|) : IKaleckil | |1945|) : iParetd l ll89'n) : lPikettv and Saed (120031 '). 
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has virtually no counterpart in modern economics, so 
econophysicists initiated a new direction here. Only 
the search theory of money ( Kivotaki and Wrightl . 
I1993I ) is somewhat related to it. This theory was 
an inspiration for the early econ ophysics paper by 
iBak. Norrelvkke. and ShubiS ( 1999t ) studying dynamics 
of money. However, a probability distribution of money 
among the agents was only recently obtained within 
the sea r ch-th eoretical approach by the economist Miguel 
iMolicd (I2006D . His distrib ution i s qu alitatively similar to 
the d i stribu tions foun d bvlAngld lfl98 6. 1992, 1993, 19jl 
I2OO2I I2OO6I ) and by llspolatov. Krapivskv. and Redned 
(|l998r ). but its functional form is unknown, because it 
was obtained only numerically. 



A. The Boltzmann-Gibbs distribution of energy 

The fundamental law of equilibrium statistical me- 
chanics is the Boltzmann-Gibbs distribution. It states 
that the probability P{e) of finding a physical system or 
subsystem in a state with the energy e is given by the 
exponential function 



Pie) 



-e/T 



(1) 



where T is the t empe rature, and c is a normalizing con- 
stant (IWannieil . [T987h . Here we set the Boltzmann con- 
stant ks to unity by choosing the energy units for mea- 
suring the physical temperature T. Then, the expecta- 
tion value of any physical variable x can be obtained as 



-Sk/T 



(2) 



where the sum is taken over all states of the system. 
Temperature is equal to the average energy per particle: 
T ^ (e), up to a numerical coefficient of the order of 1. 



Eg . IH) can be derived in different ways ( Wannierl . 
Il987l ). All derivations involve the two main ingredients: 
statistical character of the system and conservation of 
energy e. One of the shortest derivations can be sum- 
marized as follows. Let us divide the system into two 
(generally unequal) parts. Then, the total energy is the 
sum of the parts: e — ei + 62, whereas the probability 
is the product of probabilities: P{e) = P{si) P{£2)- The 
only solution of these two equations is the exponential 
function ([T]). 

A more sophisticated derivation, proposed by Boltz- 
mann, uses the concept of entropy. Let us consider N 
particles with the total energy E. Let us divide the en- 
ergy axis into small intervals (bins) of width Ae and 
count the number of particles Nk having the energies 
from Ek to Ek + Ae. The ratio N^/N = Pk gives the 
probability for a particle to have the energy Sk- Let us 
now calculate the multiplicity which is the number 
of permutations of the particles between different energy 
bins such that the occupation numbers of the bins do 



not change. This quantity is given by the combinatorial 
formula in terms of the factorials 



W = 



iV! 



(3) 



The logarithm of multiplicity is called the entropy S = 
In W. In the limit of large numbers, the entropy per 
particle can be written in the following form using the 
Stirling approximation for the factorials 



(4) 



Now we would like to find what distribution of particles 
among different energy states has the highest entropy, 
i.e., the highest multiplicity, provided the total energy of 
the system, E — NkSk, has a fixed value. Solution of 
this problem can be e asily obtained u sing the method of 
Lagrange multipHers (jWannieil . Il987l ). and the answer is 
given by the exponential distribution ([1]). 

The same result can be also derived from the er- 
godic theory, which says that the many-body system oc- 
cupies all possible states of a given total energy with 
equal probabilities . The n it is straightforward to show 
( Lopez- Ruiz et all . l2008h that the probability distribu- 
tion of the energy of an individual particle is given by 
Eq. 



B. Conservation of money 



The derivations outlined in Sec. III.AI are very general 
and only use the statistical character of the system and 
the conservation of energy. So, one may expect that the 
exponential Boltzmann-Gibbs distribution ([T]) would ap- 
ply to other statistical systems with a conserved quantity. 

The economy is a big statistical system with millions 
of participating agents, so it is a promising target for ap- 
plications of statistical mech anics. Is there a conserved 
quant ity in the economy? iDragulescu and Yakovenkol 
(|2000D argued that such a conserved quantity is money 
m. Indeed, the ordinary economic agents can only re- 
ceive money from and give money to other agents. They 
are not permitted to "manufacture" money, e.g., to print 
dollar bills. Let us consider an economic transaction be- 
tween agents i and j. When the agent i pays money 
Am to the agent j for some goods or services, the money 
balances of the agents change as follows 



rrii 
to' 



rrii 

= TO., 



Ato, 
h Am. 



(5) 

The total amount of money of the two agents before and 
after transaction remains the same 

nii + nij — ra[ -f to^- , (6) 

i.e., there is a local conservation law for money. The rule 
for the transfer of money is analogous to the trans- 
fer of energy from one molecule to another in molecular 
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collisions in a gas, and Eq. ([6]) is analogous to conser- 
vation of energy in such collisions. Conservative models 
of this kind are also studied in some econo mic literature 
(|Kivotaki and Wright! . [Tool iMolicoL l2006^ . 

We should emph asize that, in the model of 
iDragulescu and Yako venko (l200Clh [as i n the e c onom ic 
models of iKivotaki and Wrightl llQQSh : iMolicd (|2006h ]. 
the transfer of money from one agent to another repre- 
sents payment for goods an d services in a market econ- 
omy. However, the model of IDragulescu and Yakovenko! 
(j2000) only keeps track of money flow, but does not keep 
track of what goods and service are delivered. One reason 
for this is that many goods, e.g., food and other supplies, 
and most services, e.g., getting a haircut or going to a 
movie, are not tangible and disappear after consump- 
tion. Because they are not conserved, and also because 
they are measured in different physical units, it is not 
very practical to keep track of them. In contrast, money 
is measured in the same unit (within a given country 
with a single currency) and is conserved in local transac- 
tions ([6]) , so it is straightforward to keep track of money 
flow. It is also important to realize that an increase in 
material production does not produce an automatic in- 
crease in money supply. The agents can grow apples on 
trees, but cannot grow money on trees. Only a central 
ban k has the monopoly of chan ging the monetary base 
Mb (jMcConnell and Bruel . [1996!) . (De bt and credit issues 
are discussed separately in Sec. III.Dl ) 

Unlike, ordinary economic agents, a central bank or a 
central government can inject money into the economy, 
thus changing the total amount of money in the system. 
This process is analogous to an influx of energy into a sys- 
tem from external sources, e.g., the Earth receives energy 
from the Sun. Dealing with these situations, physicists 
start with an idealization of a closed system in thermal 
equilibrium and then generalize to an open system sub- 
ject to an energy flux. As long as the rate of money 
influx from central sources is slow compared with relax- 
ation processes in the economy and does not cause hy- 
perinflation, the system is in quasi-stationary statistical 
equilibrium with slowly changing parameters. This sit- 
uation is analogous to heating a kettle on a gas stove 
slowly, where the kettle has a well-defined, but slowly in- 
creasing, temperature at any moment of time. A flux of 
money may be also produced by international transfers 
across the boundaries of a country. This process involves 
complicated issues of m ultiple currencies in the world 
and their exchange rates (jMcCauie3,[lQ0i). Here we use 
an idealization of a closed economy for a single country 
with a single currency. Such an idealization is common 
in economic literature. For examp le, in the two-volume 
Hand book of Monetary Economics ( Friedman and Hahn! . 
llQQOf ) , only the last chapter out of 23 chapters deals with 
an open economy. 

Another potential problem with conservation of money 
is debt. This issue will be discussed in S ec. III.Dl As a 
starting point . [Pragulescu and Yakovenko! ( |2000) consid- 
ered simple models, where debt is not permitted, which 



is also a common idealization in some econom ic liter- 
ature (|Kivotaki and Wrightl . 11993!: iMohcol . I2006D . This 
means that money balances of the agents cannot go be- 
low zero: > for all i. Transaction ([5]) takes place 
only when an agent has enough money to pay the price: 
rrii > Am, otherwise the transaction does not take place. 
If an agent spends all money, the balance drops to zero 
rUi = 0, so the agent cannot buy any goods from other 
agents. However, this agent can still receive money from 
other agents for delivering goods or services to them. In 
real life, money balance dropping to zero is not at all 
unusual for people who live from paycheck to paycheck. 

Enforcement of the local conservation law ^ is the 
key feature for successful functioning of money. If the 
agents were permitted to "manufacture" money, they 
would be printing money and buying all goods for noth- 
ing, which would be a disaster. The physical medium of 
money is not essential here, as long as the local conser- 
vation law is enforced. The days of gold standard are 
long gone, so money today is truly the fiat money, de- 
clared to be money by the central bank. Money may 
be in the form of paper currency, but today it is more 
often represented by digits on computerized bank ac- 
counts. The local conservation law ^ is consistent with 
the fundamental principles of accounting, whether in the 
single-entry or the double-entry form. More discussion of 
banks, debt, and credit will be given in Sec. III.Dl How- 
ever, the macroeconomic monetary policy issues, such as 
mone y supply and money demand ( Friedman and Hahn! . 
|1990D . are outside of the scope of this paper. Our goal 
is to investigate the probability distribution of money 
among economic agents. For this purpose, it is appropri- 
ate to make the simplifying macroeconomic idealizations, 
as described above, in order to ensure overall stability of 
the system and existence of statistical equilibrium in the 
model. The concept of "equilibrium" is a very common 
idealization in economic literature, even though the real 
economies might never be in equilibrium. Here we ex- 
tend this concept to a statistical equilibrium, which is 
characterized by a stationary probability distribution of 
money P(m), as opposed to a mechanical equilibrium, 
where the "forces" of demand and supply match. 



C. The Boltzmann-Gibbs distribution of money 

Having recognized the principle of local money con- 
servation, IDragulescu and Yakovenko! (|200Cl[ ) argued that 
the stationary distribution of money P{m) should be 
given by the exponential Boltzmann-Gibbs function anal- 
ogous to Eq. ^ 



P{m) 



ce 



(7) 



Here c is a normalizing constant, and Tm is the "money 
temperature" , which is equal to the average amount of 
money per agent: T = (m) = M/N, where M is the total 



6 



money, and N is the number of agents.^ 

To verify this conjecture, iDragulescu and Yakovenkol 
(|2000D performed agent-based computer simulations of 
money transfers between agents. Initially all agents were 
given the same amount of money, say, $1000. Then, a 
pair of agents was randomly selected, the amount 
Am was transferred from one agent to another, and the 
process was repeated many times. Time evolution of 
the probability distribution of mon ey P(m) is illustrated 
in computer animation video s by IChen and Yakovenkol 
|2007) and by Wrigld ^00^. After a transitory pe- 
riod, money distribution converges to the stationary form 
shown in Fig. [TJ As expected, the distribution is well fit- 
ted by the exponential function ([7]). 

Several different rules f o r Am were considered by 
IDragules cu and Yakovenkol (|2000h . In one model, the 
transferred amount was fixed to a constant Am = $1. 
Economically, it means that all agents were selling their 
product s for the same price Am = $1. Computer ani- 
mation (jChen and Yakovenkol l2007t ) shows that the ini- 
tial distribution of money first broadens to a symmet- 
ric Gaussian curve, characteristic for a diffusion pro- 
cess. Then, the distribution starts to pile up around 
the TO = state, which acts as the impenetrable bound- 
ary, because of the imposed condition m > 0. As a re- 
sult. Pirn) becomes skewed (asymmetric) and eventu- 
ally reaches the stationary exponential shape, as shown 
in Fig. [1] The boundary at m = is analogous to 
the ground-state energy in statistical physics. Without 
this boundary condition, the probability distribution of 
money woul d not reach a sta tion ary st a te. Computer 
animations (|Chen and Yakovenkol l2007t ImTghtl . l2007l ) 
also show how the entropy of money distribution, defined 
as S/N = — ^^,P{mk)\'a.P{mk), grows from the initial 
value 5 = 0, where all agents have the same money, to 
the maximal value at the statistical equilibrium. 

While the model with Am = 1 is very simple and 
instructive, it is not realistic, because all prices are 
taken to be the same. In anothe r model considered by 
IDragulescu and Yakovenkol ( 2000t ). Am in each transac- 
tion is taken to be a random fraction of the average 
amount of money per agent, i.e., Ato = v{M/N), where 
J/ is a uniformly distributed random number between 
and 1. The random distribution of Ato is supposed to 
represent the wide variety of prices for different products 
in the real economy. It reflects the fact that agents buy 
and consume many different types of products, some of 
them simple and cheap, some sophisticated and expen- 
sive. Moreover, different agents like to consume these 
products in different quantities, so there is a variation in 
the paid amounts Ato, even when the unit price of the 
same product is constant. Computer simulation of this 
model produces exactly the same stationary distribution 
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* Because debt is not permitted i n this model, we have M = Mi,, 
where Mf, is the monetary base llMcConnell and Brud . Il996l) . 



FIG. 1 Histogram and points: Stationary probability dis- 
tribution of money P{m) obtained in agent-based computer 
simulations. Solid curves: Fits to the Boltzmann-Gibbs law 
(|7[). Vertical line: The ini t ial di stribution of money. From 
IDragulescu and Yakovenkol (|2000D . 



([7|), as in the first mo del. Compu ter animation for this 
model is also given by IChen andY akovcnko ( 2007 ) . 

The final distribution is universal despite differ- 
ent rules for Am. To amplify this point further, 
I Dragulescu and Yakovenkol ([2000) also considered a toy 
model, where Am was taken to be a random fraction 
of the average amount of money of the two agents: 
Am = firrii + mj)/2. This model produced the same 
stationary distribution ([7|) as the two other models. 

The models of pairwise money transfer are attrac- 
tive in their simplicity, but they represent a rather 
primitive market. Modern economy is dominated 
by big firms, which consist of many agents, so 
IDragulescu and Yakovenkol (l200(l also studied a model 
with firms. One agent at a time is appointed to be- 
come a "firm" . The firm borrows capital K from another 
agent and returns it with interest hK , hires L agents and 
pays them wages uj, manufactures Q items of a prod- 
uct, sells them to Q agents at a price p, and receives 
profit F = pQ — LuL — hK. All of these agents are ran- 
domly selected. The parameters of the model are opti- 
mized following a procedure from economics textbooks 
( McConnell and Brud . Il99d) . The aggregate demand- 
supply curve for the product is taken in the form p{Q) — 
v/Q^, where Q is the quantity consumers would buy at 
the price p, and r] and v are some parameters. The pro- 
duction function of the firm has the traditional Cobb- 
Douglas form: Q{L, K) ~ L^K^~^, where x is a param- 
eter. Then the profit of the firm F is maximized with 
respect to K and L. The net result of the firm activ- 
ity is a many-body transfer of money, which still satis- 
fies the conservation law. Computer simulation of this 
model generates the same exponential distribution ([7|), 
independently of the model parameters. The reasons for 
the universality of the Boltzmann-Gibbs distribution and 
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its limitations are discussed i n Sec. ITl.FI 

Wel l after the paper by iDragulescu and Yakovenkol 
(I2OOOI) appea r ed, the Italian econophysicists 
iPatriarca et al ] 12005) found that similar ideas had 
been pub lished ea r lier i n obsc ure Italian journals by 
Eleonora iBennatil ( 19881 Il993l ). It was proposed to 
call t hese models th e Bennati-DraRules c u-Yak ovenko 
game (|Garibaldi et all l2007t IScalas et all l2006t) . The 
Boltzmann distribution was independ ently applied t o 
social sciences by t he ph ysicist Jiirgen iMimkesI ( 2000l) : 
iMimkes and Willis! 

using the Lagrange principle 
of maximization with constraints. The exponential 
distribu tion of r noney was also found by the economist 
Martin IShubikI (|l999l ) using a Markov chain approach 
to strategic market games. A long time ago, Bcnoit 
iMandelbrotI (|l960l . p 83) observed: 



"There is a great temptation to consider the 
exchanges of money which occur in economic 
interaction as analogous to the exchanges of 
energy which occur in physical shocks be- 
tween gas molecules." 

He realized that this process should result in the expo- 
nential distribution, by analogy with the barometric dis- 
tribution of density in the atmosphere. However, he dis- 
carded this idea, because it does not produce the Pareto 
power law, and proceeded to study the stable Levy distri- 
butions. Ironically, the actual economic data, discussed 
in Sees. IIII.CI and IIV.AI do show the exponential distri- 
bution for the majority of the population. Moreover, the 
data have a finite variance, so the stable Levy distribu- 
tions are not applicable because of their infinite variance. 



D. Models with debt 

Now let us discuss how the results change when debt is 
permitted.^ From the standpoint of individual economic 
agents, debt may be considered as negative money. When 
an agent borrows money from a bank (considered here as 
a big reservoir of money), ^ the cash balance of the agent 
(positive money) increases, but the agent also acquires 
a debt obligation (negative money), so the total balance 
(net worth) of the agent remains the same. Thus, the act 
of borrowing money still satisfies a generalized conserva- 
tion law of the total money (net worth), which is now 
defined as the algebraic sum of positive (cash M) and 



^ The i deas presented here are quite similar to those by ISoddvl 
lll933t) . 

® Here we treat the bank as being outside of the system consisting 
of ordinary agents, because we are interested in money distribu- 
tion among these agents. The debt of agents is an asset for the 
bank, and deposits of cash in to the bank are liabilities of the bank 
llMcConnell and Brug.ll996D . We do not go into these details in 
order to kee p our presentation simple. For more discussion, see 
iKeenI ll2008t) . 



negative (debt D) contributions: M — D = Mj,. After 
spending some cash in binary transactions ([5]), the agent 
still has the debt obligation (negative money) , so the to- 
tal money balance of the agent (net worth) becomes 
negative. We see that the boundary condition rrii > 0, 
discussed in Sec. IILB) does not apply when debt is per- 
mitted, so m = is not the ground state any more. The 
consequence of permitting debt is not a violation of the 
conservation law (which is still preserved in the gener- 
alized form for net worth), but a modification of the 
boundary condition by permitting agents to have neg- 
ative balances mt < of net worth. A more detailed 
discussion of positive and negative money and the book- 
keeping accounting from the econoph ysics point of view 
was presented by t he physicist Dieter iBraunI ( 200ll ) and 
iFischer and BraunI (|2003a .b). 

Now we can repeat the simulation described in Sec. 
III. CI without the boundary condition m > by allowing 
agents to go into debt. When an agent needs to buy a 
product at a price Am exceeding his money balance m^, 
the agent is now permitted to borrow the difference from 
a bank and, thus, to buy the product. As a result of 
this transaction, the new balance of the agent becomes 
negative: = — Am < 0. Notice that the local con- 
servation law ([5]) and ([6]) is still satisfied, but it involves 
negative values of m. If the simulation is continued fur- 
ther without any restrictions on the debt of the agents, 
the probability distribution of money P{rn) never stabi- 
lizes, and the system never reaches a stationary state. As 
time goes on, P{m) keeps spreading in a Gaussian man- 
ner unlimitedly toward m = +00 and m — —00. Because 
of the generalized conservation law discussed above, the 
first moment (m) = M\j/N of the algebraically defined 
money m remains constant. It means that some agents 
become richer with positive balances m > at the ex- 
pense of other agents going further into debt with nega- 
tive balances to < 0, so that M — Mb H- D. 

Common sense, as well as the experience with the cur- 
rent financial crisis, tells us that an economic system can- 
not be stable if unlimited debt is permitted.'' In this case, 
agents can buy any goods without producing anything in 
exchange by simply going into unlimited debt. Arguably, 
the current financial crisis was caused by the enormous 
debt accumulation in the system, triggered by subprime 
mortgages and financial derivatives based on them. A 
widely expressed opinion is that the current crisis is not 
the problem of liquidity, i.e., a temporary difficulty in 
cash flow, but the problem of insolvency, i.e., the inher- 
ent inability of many participants pay back their debts. 

Detailed discussion of the current economic situation 
is not a subject of this paper. Going back to the idealized 
model of money transfers, one would need to impose some 
sort of modified boundary conditions in order to prevent 



In qu alitatively agreement with the conclusions by iMcCauIevI 
l l2008t) . 
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N=500, M=5*10 , time=4*10^. 

I8n 1 ' ' ^ 
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Money, m 

FIG. 2 Histograms: Stationary distributions of money with 
and without debt. The debt is limited to rud = 800. 
Solid curves: Fits to the Boltzmann-Gibbs laws with the 
"money temperatures" Tm = 18 00 and T,„ — 1000. From 
iDraeulescu and Yakovenkol (|2000l ) . 



150 

Monetary Wealth, m 

FIG. 3 The stationary distribution of money for the required 
reserve ratio R = 0.8. The distribution is exponential for 
positive and negative money with different "temperatures" 
r+ a nd T-, as illustrated by t he inset on log-linear scale. 
From IXi. Ding, and Wand (|2005l ) . 



unlimited growt h of debt and to ensure ove r all sta bility 
of the system. iDragulescu and Yakovenkol (2Q03) con- 
sidered a simple model where the maximal debt of each 
agent is limited to a certain amount nid- This means 
that the boundary condition mi > is now replaced by 
the condition > —nid for all agents i. Setting inter- 
est rates on borrowed money to be zero for simplicity, 
IDragulescu and Yakovenkol ( 20001 ) performed computer 
simulations of the models described in Sec. III. Cl with the 
new boundary condition. The results are shown in Fig. 
[2l Not surprisingly, the stationary money distribution 
again has the exponential shape, but now with the new 
boundary condition at m = — and the higher money 
temperature Td = rud + M^/N. By allowing agents to go 
into debt up to md, we effectively increase the amount of 
money available to each agent by rud- So, the money tem- 
perature, which is equal to the average amount of effec- 
tively available money per agent, increases correspond- 
ingly^ 

IXi. Ding, and Wand (|2005D considered another, more 
realistic boundary condition, where a constraint is im- 
posed not on the individual debt of each agent, but 
on the total debt of all agents in the system. This is 
accomplished via the req uired reserve ratio R, whic h 
is briefly explained below I McConnell and Brud . Il996t ). 
Banks are required by law to set aside a fraction R of 
the money deposited into bank accounts, whereas the 
remaining fraction 1 — R can be loaned further. If the 
initial amount of money in the system (the money base) 
is Mb, then, with repeated loans and borrowing, the to- 
tal amount of positive money available to the agents in- 
creases to Af = Mh/R, where the factor is called 

vt- 



the money multiplier ([McConnell and" Bruelll996l ). This 
is how "banks create money". Where does this extra 



money come from? It comes from the increase in the to- 
tal debt in the system. The maximal total debt is given 
hy D — Mb/ R~ Mb and is limited by the factor R. When 
the debt is maximal, the total amounts of positive, Mb/ R, 
and negative, — R)/R, money circulate among the 

agents in the system, so there are two constraints in the 
model considered bv lXi. Ding, and Wa ng ( ,2005j . Thus, 
we expect to see the exponential distributions of positive 
and negative money characterized by two different tem- 
peratures: T+ = Mb/RN and T. = Mb{l - R)/RN. 
This is exactly what was found in computer simula- 
tions by IXi. Ding, and Wand (|2005[ ). as shown in Fig. 
51 Similar two-sided dist ributions were also found by 



Fischer and BraunI (|2003al ) 



However, in reality, the reserve requirement is not 
effective in stabilizing total debt in the system, be- 
cause it applies only to dep osits from gene ral public, 
but not from corporations ( O'Brienl . l2007t ).^ More- 
over, there are alternative instruments of debt, includ- 
ing derivatives and various unregulated "financial inno- 
vations". As a result, the total debt is not limited in 
practice and sometimes can reach catastrophic propor- 
tions. Here we briefly discuss several models with non- 
stationary de bt. Thus far, we did not conside r the in- 
terest rates. IDragulescu and Yakovenkol ( 2000l ) studied 
a simple model with different interest rates for deposits 
into and loans from a bank. Computer simulations found 
that money distribution among the agents is still expo- 
nential, but the money temperature slowly changes in 



* Australia does not have reserve requirements, but China actively 
uses reserve requirements as a tool of monetary policy. 
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time. Depending on the choice of parameters, the total 
amount of money in circulation either increases or de- 
creases in time. A more sophisticated m acroeconomic 
model was studied by the economist Steve iKeenI ( 19951 
He found that one of the regimes is the debt- 
induced breakdown, where all economic activity stops 
under the burden of heavy debt and cannot be restarted 
without a "debt moratorium" . The interest rates were 
fixed in these models a n d not adjusted self-consistently. 
ICockshott and Cottrelll ( 20081 ) proposed a mechanism, 
where the interest rates are set to cover probabilistic 
withdrawals of deposits f rom a bank. In an agen t -based 
simulation of the model, ICockshott and Cottrelll ( 20081 ) 
found that money supply first increases up to a certain 
limit, and then the economy experiences a spectacular 
crash under the weight of accumulated debt. Further 
studies along these lines would be very interesting. In 
the rest of the paper, we review various models without 
debt proposed in literature. 



E. Proportional money transfers and saving propensity 



N=500, M=5*10=, a=1/3. 



In the models of money transfer discussed in Sec. Ill.Ci 
the transferred amount Am is typically independent of 
the money balances of the agents involved. A differ- 
ent model was introduced in physics liter ature earlier by 
llspolatov. Krapivskv. and Redned ( 19981 ) and called the 
multiplicative asset exchange model. This model also sat- 
isfies the conservation law, but the transferred amount of 
money is a fixed fraction 7 of the payer's money in Eq. 



Am = "frrii. 



(8) 



The stationary distribution of money in this model, com- 
pared in Fig. m with an exponential function, is similar, 
but not exactly equal, to the Gamma distribution: 



P{m) 



cra^ e 



-m/T 



(9) 



Eq. ^ differs from Eq. ([7]) by the power-law prefactor 
From the Boltzmann kinetic equation (d i scusse d in 



i/3 



Sec. HLF]) . ils polatov. Krapivskv. and Redned (|l998l ) de- 
rived a formula relating the parameters 7 and j3 in Eqs. 
(IHD and m-. 



(3^ -1 -ln2/ln(l -7). 



(10) 



When payers spend a relatively small fraction of their 
money 7 < 1/2, Eq. (fTU|) gives /? > 0. In this case, 
the population with low money balances is reduced, and 
P(0) = 0, as shown in Fig IH 

The economist Thomas i Luxl (l2005h brought to the at- 
tention of physicists that essentially the same model, 
called the inequality process, had been introduc ed and 
studied muc h earl i er by the s ociologist John lAngld 
il986, 1992. ,19931 Il996l l2002t) . see also the review 
bv lAngld (120061) for additional refere nces. While 
llspolatov. Krapivskv. and Redned ( 19981 ) did not give 
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FIG. 4 Histogram: Stationary probability distribution of 
money in the multiplicative random exchange model ^ for 
7 = 1 /3. Solid c urve: The exponential Boltzmann-Gibbs law. 



1/3 
m|D 



From iDragulescu and Yakovenkg (|2000,) . 



much justification for the proportionality law (|5|) . lAngld 
( 19861 ) connected this rule with the surplus theory of so- 
cial stratification ( Engeli . I1972D . which argues that in- 
equality in human society develops when people can pro- 
duce more than necessary for minimal subsistence. This 
additional wealth (surplus) can be transferred from origi- 
nal producers to othe r peop l e, thu s generating inequality. 
In the first paper bv I Anglel (|1986( ). the parameter 7 was 
randomly distributed, and another parameter 5 gave a 
higher probability of winning to the agent with the higher 
money balance in Eq. ([5]). However, in the following pa- 
pers, he simplified the model to a fixed 7 (denoted as uj by 
Angle) and equal probabilities of winning for higher- and 
lower-balance agents, which makes it completely equiv- 
alent t o the model of Ispolatov. Krapivskv. and Redned 
(I1998D . lAlgl3(|2002ll2006l) also considered a model where 
groups of agents have different values of 7, simulating 
the effect of education and other "human capital". All 
of these models generate a Gamma-like distribution, well 
approximated by Eq. 

Another m odel with an element of propor t ionali ty was 
proposed bv IChakraborti and Chakrabartil ( 2000( ).^ In 
this model, the agents set aside (save) some fraction of 
their money Am^ , whereas the rest of their money balance 
(1 — A)mi becomes available for random exchanges. Thus, 
the rule of exchange ([5]) becomes 

m- = XttLi + ^(1 - X){mi + urij), 
m',j ^ Xnij + (1 - 0(1 - ^)(f"'» + mj). (11) 
Here the coefficient A is called the saving propensity, and 



® This paper originally appeared as a follow-up e-print 
cond-mat/0004256 on the e-p rint |cond-mat/0001432| by 



[PriSigulescu and Yakovenkd 1 I2OOOD . 
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the random variable ^ is unifor mly dist r ibuted between 
and 1. It was pointed out by lAngl j ( 20061 ) that, by 
the change of notation A (1 — 7), Eq. ([TT|) can be 
transformed to the same form as Eq. ([5]), if the random 
variable ^ takes only discrete values and 1. Compu ter 
simulations by IChakraborti and Chakrabartil ( 2000f ) of 
the model (1111) found a stationary distribution close to 
the Gamma distribution ([5]). It was shown that the pa- 
rameter P is related t o the saving propensity A by the for- 
mula / 3 = SA/fl- A) dPat riarca. Chakraborti. and KaskiL 
2004allbl: iPatriarca et al\. 12005 : 



Repetowicz. Hutzler. and Richmondl . l2005l ). For A ^ 



agents always keep some money, so their balances never 
drop to zero, and P{0) = 0, whereas for A = the 
distribution becomes exponential. 



In the sub s equen t papers by the Kolkata school 
(IChakrabartiL I2005D and related papers, the case 
of random saving propensity was studied. In these 
models, the agents are assigned random parameters 
A drawn from a uniform distribution b e tween and 
1 (jChatteriee. Chakrabarti. and Manual . \2004 ). It 
was found that this model produces a power-law tail 
P{m) oc 1/m? at high m. The reasons for stability of 
this law wer e understood using the Boltzmann kinetic 
equat ion (jChatt eriee. Chakrabar ti. and Stinchcombe 



Das and Yarlagaddai 



20051: 

Repetowicz. Hutzler. and Richmondl 
but most e legantly i n the mean-field 
(iBhattacharvva . Chattcri ee. and Cha kraba rti 
IChatteriee and Chakrabarti 2007t iMohantvl . 
The fat tail originates from the agents whose saving 
propensity is clo se to 1, who hoard money and do not 
give i t back ( Patriarca. Chakraborti. and Germand . 




l2006t IPatriarca etaH . 
mathematical 



ous 



2OO5I ). A more rigor- 

treatment of the pro ble m wa s 

given by Diiring. Matthes. and Toscani^ ('2008'); 
iDtiring and Toscanil (2007); Mat thes and T oscani 
( 20081 ). An interesting matrix formu lation of the 
problem was presented by ICuptal ( 20061 ). Relaxation 
rate in the money transfer rnodels \ yas stu d ied b y 
Diiring. Matthes. and Toscanil (|2008|): ICuptal (l2008l): 



Patriarca et al\ (|2007l ). iDragulescu and Yakovenkol 
(|200flh considered a model with taxation, which also has 
an element of proportionality. The Gamma distribution 
was also studied for conservati ve mod el s wit hin a 
simple Boltzmann approach by iFerrerol ( 2004 ) and, 
using more complicated ru les of exchange motivated 
by pol itical economy, by IScafetta . Picozzi, an d We st 
2004al lbl). Independently, the economist Miguel iMolicol 



20061 ) studied conservative exchange models where 
agents bargain over prices in their transactions. He 
found stationary Gamma-like distributions of money in 
numerical simulations of these models. 



F. Additive versus multiplicative models 

The stationary distribution of money ([9]) for the mod- 
els of Sec. III.EI is different from the simple exponential 
formula ([7|) found for the models of Sec. III. CI The origin 
of this difference can be understood from the Boltzmann 
kinetic equation ("Lifshit z and Pitaevskiil . Il98ll : IWannierl . 
Il987l ) . This equation describes time evolution of the dis- 
tribution function P(m) due to pairwise interactions: 



dP{m) 
Jt 



m.m' |— >-|m — A,m 



_A]P(m)P(m') (12) 



+f[,n-Ajn'+A]^[m,ni']P{m - ^)Pim' + A)} dm' dA. 

Here /[m,m']^[m-A,m'-i-A] IS the probability of transfer- 
ring money A from an agent with money m to an agent 
with money m' per unit time. This probability, multi- 
plied by the occupation numbers P{m) and P{m'), gives 
the rate of transitions from the state [m, m'] to the state 
[m — A, m' -t- A]. The first term in Eq. (fT^ gives the de- 
population rate of the state m. The second term in Eq. 
(fT2| describes the reversed process, where the occupation 
number P{m) increases. When the two terms are equal, 
the direct and reversed transitions cancel each other sta- 
tistically, and the probability distribution is stationary: 
dP{m)/dt = 0. This is the principle of detailed balance. 

In physics, the fundamental microscopic equations of 
motion obey the time-reversal symmetry. This means 
that the probabilities of the direct and reversed processes 
are exactly equal: 



/[n 



-A,i 



-A] — /[m-A,m'+A]^[r7 



(13) 



When Eq. (fT5|) is satisfied, the detailed balance condi- 
tion for Eq. (fT^ reduces to the equation P{m)P{m') = 
P{m — A)P{m' + A), because the factors / cancels out. 
The only solution of this equation is the exponential func- 
tion P{m) = cexp{—ni/Tm), so the Boltzmann-Gibbs 
distribution is the stationary solution of the Boltzmann 
kinetic equation (fT^ . Notice that the transition prob- 
abilities (fT5|) are determined by the dynamical rules of 
the model, but the equilibrium Boltzmann-Gibbs distri- 
bution does not depend on the dynamical rules at all. 
This is the origin of the universality of the Boltzmann- 
Gibbs distribution. We see that it is possible to find the 
stationary distribution without knowing details of the dy- 
namical rules (which arc rarely known very well) , as long 
as the symmetry condition (|13p is satisfied. 

The models considered in Sec. III.CI have the time- 
reversal symmetry. The model with the fixed money 
transfer A has equal probabilities p3p of transferring 
money from an agent with the balance m to an agent 



with the balance 



and vice versa. This is also true 



when A is random, as long as the probability distribution 
of A is independent of m and m'. Thus, the stationary 
distribution P[m) is always exponential in these models. 

However, there is no fundamental reason to expect 
the time-reversal symmetry in economics, where Eq. (jl3p 
may be not valid. In this case, the system may have a 
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non-exponential stationary distribution or no stationary 
distribution at all. In the model ([5]), the time- reversal 
symmetry is broken. Indeed, when an agent i gives a 
fixed fraction 7 of his money rrii to an agent with bal- 
ance TTij, their balances become (1 — "f)mi and rrij + jmi. 
If we try to reverse this process and appoint the agent j 
to be the payer and to give the fraction 7 of her money, 
7(TOj -|- jrrii), to the agent i, the system does not return 
to the origin al configuration [nii, nij]. As emphasized by 
lAngld (|200d) . the payer pays a deterministic fraction of 
his money, but the receiver receives a random amount 
from a random agent, so their roles are not interchange- 
able. Because the proportional rule typically violates 
the time-reversal symmetry, the stationary distribution 
P{m) in multiplicative models is typically not exponen- 
tial. Making the transfer dependent on the money bal- 
ance of the payer effectively introduces Maxwell's demon 
into the model. Another view on the time-reversal sym- 
metry in economic dynamics was presented by Ao (2007|. 

These examples show that the Boltzmann-Gibbs dis- 
tribution does not necessarily hold for any conservative 
model. However, it is universal in a limited sense. For 
a broad class of models that have time-reversal symme- 
try, the stationary distribution is exponential and does 
not depend on details of a model. Conversely, when 
the time-reversal symmetry is broken, the distribution 
may depend on details of a model. The difference be- 
tween these two classes of models may be rather sub- 
tle. Deviations from the Boltzmann-Gibbs law may oc- 
cur only if the transition rates / in Eq. (fTS]) explicitly 
depend on the agents' money m or m' in a n asymmetric 
manner. iDragulescu and Yakovenkol ( 2000l ) performed a 
computer simulation where the direction of payment was 
randomly fixed in advance for every pair of agents (i, j). 
In this case, money flows along directed links between 
the agents: « — > j — > fc, and the time-reversal symme- 
try is strongly violated. This model is closer to the real 
economy, where one typically receives money from an 
employer and pays it to a grocery store. Nevertheless, 
the Boltzmann-Gibbs distribution was still found in this 
model, because the transition rates / do not explicitly de- 
pend on m and to' and do not violate Eq. (|13p . A more 
general study of money exchange m o dels o n directed net- 
works was presented bv lChatteried ( 20091 ). 



In the absence of detailed knowledge of real mi- 
croscopic dynamics of economic exchanges, the semi- 
universal Boltzmann-Gibbs distribution ([7]) is a nat- 
ural starting point. Moreover , the assumption of 
IDragulescu and Yakovenkol ( 2000( ) that agents pay the 
same prices Ato for the same products, independent of 
their money balances m, seems very appropriate for the 
modern anonymous economy, especially for purchases 



over the Internet. There is no particular empirical ev- 
idence for the proportional rules ^ or pT|) . However, 
the difference between the additive ([7]) and multiplicative 
© distributions may be not so crucial after all. From the 
mathematical point of view, the difference is in the im- 
plementation of the boundary condition at to = 0. In 
the additive models of Sec. III. CI there is a sharp cutoff 
for P{m) 7^ at TO = 0. In the multiplicative models of 
Sec. III. El the balance of an agent never reaches to = 0, so 
P(to) vanishes at to —> in a power-law manner. But for 
large to, P(to) decreases exponentially in both models. 

By further modifying the rules of money transfer and 
introducing more parameters in the mod els, one can ob- 
tain even more complicat e d dist ributions ( Saif and Gadel . 
120071: IScafetta and Westl . I2007D . However, one can ar- 
gue that parsimony is the virtue of a good mathematical 
model, not the abundance of additional assumptions and 
parameters, whose correspondence to reality is hard to 
verify. 



III. STATISTICAL MECHANICS OF WEALTH 
DISTRIBUTION 

In the econophysics literature on exchange models, the 
terms "money" and "wealth" are often used interchange- 
ably. However, economists emphasize the difference be- 
tween these two concepts. In this section, we review the 
models of wealth distribution, as opposed to money dis- 
tribution. 



A. Models with a conserved commodity 

What is the difference betwee n money and wealth? 
IDragulescu and Yakovenkol ( 2000( ) argued that wealth Wi 
is equal to money rrii plus the other property that an 
agent i has. The latter may include durable material 
property, such as houses and cars, and financial instru- 
ments, such as stocks, bonds, and options. Money (paper 
cash, bank accounts) is generally liquid and countable. 
However, the other property is not immediately liquid 
and has to be sold first (converted into money) to be 
used for other purchases. In order to estimate the mon- 
etary value of property, one needs to know its price p. 
In the simplest model, let us consider just one type of 
property, say, stocks s. Then the wealth of an agent i is 
given by 



rrii +p Si 



(14) 



However, when Am is a fraction of the total money mi+rrij of the 
two agents, the model is time-reversible and has the exponential 
distribution, as discussed in Sec. III. CI 



It is assumed that the price p is common for all agents 
and is established by some kind of market process, such 
as an auction, and may change in time. 

It is reasonable to start with a model where both 
the total money M = rrii and the to tal stock 
S = yi. Si are conserved (lAusloos and Pckalski . 



20071: IChakraborti. Pradhan. and Chakrabartil . 12001 



Chatteriee and Chakrabartl " 2006t )" The agents pay 
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money to buy stock and sell stock to get money, 
and so on. Although M and S are conserved, the 
total wealth W = yi,- Wj i s gene rally not conserved 



( Chatteriee and Chakrabartil l2006l ). because of price 
fluctuation in Eq. p4|) . This is an important difference 
from the money transfers models of Sec. |TT1 The wealth 
Wi of an agent i, not participating in any transactions, 
may change when transactions between other agents 
establish a new price p. Moreover, the wealth Wi of an 
agent i does not change after a transaction with an agent 
j. Indeed, in exchange for paying money Am, the agent 
i receives the stock As — Am/p, so her total wealth 
(fT4)) remains the same. Theoretically, the agent can 
instantaneously sell the stock back at the same price and 
recover the money paid. If the price p never changes, 
then the wealth Wi of each agent remains constant, 
despite transfers of money and stock between agents. 

We see that redistribution of wealth in this model 
is directly related to price fluctuations. A math- 
ematical model of this proce ss was studied by 
ISilver. Slud. and Takamotd (|2002l) . In this model, the 
agents randomly change preferences for the fraction of 
their wealth invested in stocks. As a result, some agents 
offer stock for sale and some want to buy it. The price p 
is determined from th e market-clearing auction matchini 
supply and demand. ISilver. Slud. and Takamotol ( 200 
demonstrated in computer simulations and proved ana- 
lytically using the theory of Markov processes that the 
stationary distribution P{'w) of wealth w in this model 
is given by the Gamma distribution, as in Eq. 
Various modifications of this model considered by ILuxI 
[2005,1 . such as introducing monopolistic coalitions, do 
not change this result significantly, which shows robust- 
ness of the Gamma distribution. For models with a con - 
served commodity, IChatteriee and Chakrabartil ( 20061 ) 
found the Gamma distribution for a fixed saving propen- 
sity and a power-law tail for a distributed saving propen- 
sity. 

Another model with co n served money and stock was 
studied by iRaberto et al. (2003) for an artificial stock 
market, where traders follow different investment strate- 
gies: random, momentum, contrarian, and fundamen- 
talist. Wealth distribution in the model with random 
traders was found have a power-law tail P{w) ^ l/vo^ 
for large w. However, unlike in other simulations, where 
all agents initially have equal balances, here the initial 
money and stock balances of the agents were randomly 
populated according to a power law with the same ex- 
ponent. This raises the question whether the observed 
power-law distribution of wealth is an artifact of the ini- 
tial conditions, because equilibration of the upper tail 
may take a very long simulation time. 



B. Models with stochastic growth of wealth 

Although the total wealth W is not exactly con- 
served in the models considered in Sec. IIII.Al never- 



theless W remains constant on average, because the 
total money M and stock S are conserved. A dif- 
ferent model for wealth dist ribution was proposed by 
iBouchaud and MezardI ( 2000l ). In this model, time evo- 
lution of the wealth Wi of an agent i is given by the 
stochastic differential equation 



dwi . , 



mi) 



(15) 



where r]i {t) is a Gaussian random variable with the mean 
(77) and the variance 2a^. This variable represents growth 
or loss of wealth of an agent due to investment in stock 
market. The last two terms describe transfer of wealth 
between different agents, which is taken to be propor- 
tional to the wealth of the payers with the coefficients . 
So, the model ([T5|) is multiplicative and invariant under 
the scale transformation Wi — > Zwi. For simplicity, the 
exchange fractions are taken to be the same for all agents: 
Jij — J/N for all i ^ j, where N is the total number of 
agents. In this case, the last two terms in Eq. (fT5)) can 
be written as J{{w) — Wi), where {w) = '^iWi/N is the 
average wealth per agent. This case represents a "mean- 
field" model, where all agents feel the same environment. 
It can be easily shown that the average wealth increases 
in time as {w)t — (w}oe^^''^^'^ Then, it makes more 
sense to consider the relative wealth Wi = Wi/{w)t- Eq. 
(fist for this variable becomes 



dwi 

~dr 



= - {ri) - CT^) Wi + J(l - Wi). 



(16) 



The probability distribution P{w, t) for the stochastic 
differential equation (|16p is governed by the Fokker- 
Planck equation 

dP _ d[J{w-l)^'(j'^w]P ^ ^2 d fj. d{wP) \ 
dt dw dw \ dw J 

The stationary solution [dP/dt = 0) of this equation is 
given by the following formula 



P{w) 



(18) 



The distribution (fT8| is quite different from the 
Boltzmann-Gibbs ([7]) and Gamma ^ distributions. Eq. 
(fT8)) has a power- law tail at large w and a sharp cutoff at 
small w. Eq. (fT5)) is a version of the generalized Lotka- 
Volterra model, a nd the stationary distribution ( T8l) was 

2002) . 



also obtained by ISolomon and Richmond! ( 2001 



The model w as generalized to include negative wealth by 

iHuane (2004*). 

IBouchaud a nd MezardI ( 2000h used the mean-field ap- 
proach. A similar result was found fo r a model wit h 
pairwise interaction between agents by ISlanin 3 (l2004l) . 
In his model, wealth is transferred between the agents 
following the proportional rule ([5]), but, in addition, the 
wealth of the agents increases by the factor 1 -t- C in each 
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transaction. This factor is supposed to reflect creation 
of wealth in economic interactions. Because the total 
wealth in the system increases, it makes sense to con- 
sider the distribution of rela tive wea l th P( w). In the 
limit of continuous tradina;, ISlanin i (|2004l ) found the 
same stationary distribution (fT8|). This result was re- 
produced using a mathe matically more involved treat- 
ment o f this model by Cordier. Pareschi. and Toscanil 
( 20051 ): iPareschi and Toscanil (|2006 ). Numerical sim- 
ulations of the models with stochastic noise 77 by 
IScafetta. Picozzi. and West! ((2004a^b) also found a power 
law tail for large w. Equivalence between the models with 
pairwise e xchange and exchange with a reservoir was dis- 
cussed bv lBasu and Mohantvl (|2008h . 

We now contrast the models discussed in Sees. IIII.AI 
and IIII.BI In the former case, where money and com- 
modity are conserved, and wealth does not grow, the 
distribution of wealth is given by the Gamma distribu- 
tion with the exponential tail for large w. In the latter 
models, wealth grows in time exponentially, and the dis- 
tribution of relative wealth has a power-law tail for large 
w. These results suggest that the presence of a power- 
law tail is a nonequilibrium effect that requires constant 
growth or inflation of the economy, but disappears for a 
closed system with conservation laws. 

The discussed models were reviewed 
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the popular article by 
of lack of space, we 
w ith w ealth condensa- 
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Braun. 2006: 



Burda et al. . 20021: llspolatov. Krapivskv. and Redner 



19981: iPianegonda et aZ.l . l2003l ). where a few agents accu- 
mulate a finite fraction of the total wealth, and studies o f 
wealth distribution on complex networks (ICoelhq et al 



20051; iDi Matteo. Aste. and Ha^ . |2004 IHu et all 12006 , 
2007; Ig lesias et all 2003 ). So far, we discussed the 



models with long-range interaction, where any agent can 
exchange money and wealth with any other agent. A 
local model, where agents trade only with the nearest 
neighb ors, was studied by iBak. N0rrelvkke. and ShubikI 
(|1999D . 



C. Empirical data on money and wealth distributions 

It would be interesting to compare theoretical results 
for money and wealth distributions in various models 
with empirical data. Unfortunately, such empirical data 
are difficult to find. Unlike income, which is discussed in 
Sec. llVi wealth is not routinely reported by the majority 
of individuals to the government. However, in some coun- 
tries, when a person dies, all assets must be reported for 
the purpose of inheritance tax. So, in principle, there ex- 
ist good statistics of wealth distribution among dead peo- 
ple, which, of course, is different from the wealth distri- 
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FIG. 5 Cumulative probability distribution of net wealth in 
the UK shown on log-log (main panel) and log-linear (inset) 
scales. Points represent the data from the Inland Revenue, 
and solid lines are fits to the exp onential (Boltzmann-Gibbs) 
and po wer (Pareto) laws. From iDraeulescu and Yakovenkd 
(|2001bl ). 



bution among the living. Using an adjustment procedure 
based on the age, gender, and other characteristics of the 
deceased, the UK tax agency, the Inland Revenue, recon- 
structed the wealth distribution of the whole populatio n 



structea tne weaitn distribution 01 tne wnole populatio n 
of the UK ( Her Majesty Revenue and Customa . l2003f ). 



Fig. El shows the UK dat a for 19 96 reproduced from 
iDragulescu and Yakovenkd (l20fllbl ). The figure shows 
the cumulative probability C(u') ~ P(u'') dw' as a 
function of the personal net wealth w, which is composed 
of assets (cash, stocks, property, household goods, etc.) 
and liabilities (mortgages and other debts) . Because sta- 
tistical data are usually reported at non-uniform intervals 
of w, it is more practical to plot the cumulative proba- 
bility distribution C{w) rather than its derivative, the 
probability density P{w). Fortunately, when P{w) is an 
exponential or a power-law function, then C{w) is also 
an exponential or a power-law function. 

The main panel in Fig. [5] shows a plot of C{w) on the 
log-log scale, where a straight line represents a power- 
law dependence. The figure shows that the distribu- 
tion follows a power law C{w) oc 1/w" with the expo- 
nent a = 1.9 for the wealth greater than about 100 k£. 
The inset in Fig. [5] shows the same data on the log- 
linear scale, where a straight line represents an expo- 
nential dependence. We observe that, below 100 k£, 
the data are well fitted by the exponential distribution 
C{w) oc exp(— lu/Tu,) with the effective "wealth temper- 
ature" = 60 k£ (which corresponds to the median 
wealth of 41 'k£). So, the distribution of wealth is char- 
acterized by the Pareto power law in the upper tail of 
the distribution and the exponential Boltzmann-Gibbs 
law in the lower part of the distribution for the great 
majority (about 90%) of the population. Similar results 
are found for the distribution of income, as discussed in 
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Sec. IIVI One may speculate that wealth distribution in 
the lower part is dominated by distribution of money, be- 
cause the co rresponding; people do not have other signif- 
icant assets ( Levy and Levvl . l2003fl . so the results of Sec. 
im give the Boltzmann-Gibbs law. On the other hand, 
the upper tail o f wealth distribution is dominated by in- 
vestment assess ( Lew and Levvl [20031 ) . where the results 
of Sec. IIII.BI give the Pareto l aw. The power law was 
studied by many researchers (iKlass et all . l2007t iLevvl . 
120031: iLevv and Levvl . 120031 : ISinhal . I2006D for the upper- 
tail data, such as the Forbes list of 400 richest people. On 
the other hand, statistical surveys o f the population, such 
as th e Survey of Consumer Finance ( Diaz-Gimenez et all , 
Il997t ) and the Panel Study of Income Dynamics (PSID), 
give more information ab out the lowe r part o f the wealth 
distribution. Curiously, lAbul-Magdl (|2002l ) found that 
the wealth distrib ution in the ancient Egypt was consis- 
tent with Eq. (HHI). iHegvi et all (|2007t ) found a power-law 
tail for the wealth distribution of aristocratic families in 
medieval Hungary. 

For direct comparison with the results of Sec. HIl it 
would be interesting to find data on the distribution of 
money, as opposed to the distribution of wealth. Making 
a reasonable assumption that most people keep most of 
their money in banks, one can approximate the distribu- 
tion of money by the distribution of balances on bank 
accounts. (Balances on all types of bank accounts, such 
as checking, saving, and money manager, associated with 
the same person should be added up.) Despite imperfec- 
tions (people may have accounts in different banks or 
not keep all their money in banks), the distribution of 
balances on bank accounts would give valuable informa- 
tion about the distribution of money. The data for a 
large enough bank would be representative of the distri- 
bution in the whole economy. Unfortunately, it has not 
been possible to obtain such data thus far, even though 
it would be completely anonymous and not compromise 
privacy of bank clients. 

The data on the distribution of bank accounts balances 
would be useful, e.g., to the Federal Deposits Insurance 
Company (FDIC) of the USA. This government agency 
insures bank deposits of customers up to a certain max- 
imal balance. In order to estimate its exposure and the 
change in exposure due to a possible increase in the limit, 
FDIC would need to know the probability distribution of 
balances on bank accounts. It is quite possible that FDIC 
may already have such data. 

Measuring the probability distribution of money would 
be also very useful for determining how much people 
can, in principle, spend on purchases (without going into 
debt). This is different from the distribution of wealth, 
where the property component, such as a house, a car, 
or retirement investment, is effectively locked up and, in 
most cases, is not easily available for consumer spending. 
Thus, although wealth distribution may reflect the dis- 
tribution of economic power, the distribution of money 
is more relevant for immediate consumption. 
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FIG. 6 Cumulative probability distribution of tax returns for 
USA in 1997 shown on log-log (main panel) and log-linear 
(inset) scales. Points represent the Internal Revenue Service 
data, and solid l ines are fits to the exponenti al and power-law 
functions. From lOragulescu and Yakovenkol (20031. 



IV. DATA AND MODELS FOR INCOME DISTRIBUTION 

In contrast to money and wealth distributions, more 
empirical data are available for the distribution of in- 
come r from tax agencies and population surveys. In 
this section, we first present empirical data on income 
distribution and then discuss theoretical models. 



A. Empirical data on income distribution 

Empirical studies of income distribution have a long 
history in the economic literature. ^"'^ Many articles 
on this subject appear in the journal Review of In- 
come and Wealth, published on behalf of the In- 
ternational Association for Rese arch in Incom e and 
Wealth. Following the work by iParetol ( 18971 ). much 
attention was focused on the power-law upper tail 
of income distribution and less on the lower part. 
In contrast to more com plicated functions discussed 
in the economic literature (lAtkinson and Bourguignonl 



'2000: Cha mpernowne and Cowellll 19981: lKakwanilll980l) . 
iDragulescu and Yakovenkol (|2001al ) demonstrated that 



the lower part of income distribution can be well 
fitted with the simple exponential function P{r) — 
cexp(— r/Tr), which is characterized by just one 
parameter, the "income temper a ture" T^.. Then 
IDragulescu and Yak ovenko' (l2001bl . |2003| ) showed that 
the whole income distribution can be fitted by an expo- 
nential function in the lower part and a power-law func- 



See, 



lAtkinson and Bourguignon] ||2000D: 



lAtkinson a nd Pikcttv (2007); Champer nowne and Cowelll 
l ll998t ): lKakwanii (ilQSQ ): iPikettv and Saea (.200^ 
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US, IRS data for 1983 and 2000 
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FIG. 7 Cumulative probability distribution of tax returns 
plotted on log-log scale versus r/Tr (the annual income r nor- 
malized by the average income T,- in the exponential part of 
the distribution). The IRS data points are for 1983-2001, and 
the columns of numb ers give the values of Tr for the corre- 
sponding years. From lSilva and Yakovenkol (|2005l ). 



tion in the upper part, as shown in Fig. [51 The straight 
line on the log-linear scale in the inset of Fig. [6] demon- 
strates the exponential Boltzmann-Gibbs law, and the 
straight line on the log-log scale in the main panel illus- 
trates the Pareto power law. The fact that income distri- 
bution consists of two distinct part s reveals the two-class 
structure of the American society ( Silva and Yakovenkol 
120051 : lYakovenko and Silval . 12005) . Coexistence of the ex- 
ponential and power-law distributions is also known in 
plasma physics and astrophysics, where t hey are called 
the "therm al" and "supert hermal " part s ( CoUieil . 120041 : 
iDesai et~al... . 200'i:. Hasega wa et aLl . ll985l ). The boundary 
between the lower and upper classes can be defined as the 
intersection point of the exponential and power-law fits 
in Fig.[6l For 1997, the annual income separating the two 
classes was about 120 k$. About 3% of the population 
belonged to the upper class, and 97% belonged to the 
lower class. 

ISilya and Yakovenkol (I2005D studied time evolution of 
income distribution in the USA during 1983-2001 using 
the data from the Internal Revenue Service (IRS), the 
government tax agency. The structure of income dis- 
tribution was found to be qualitatively the same for all 
years, as shown in Fig. [T] The average income in nom- 
inal dollars has approximately doubled during this time 
interval. So, the horizontal axis in Fig. [7] shows the nor- 
malized income r/T^, where the "income temperature" 
Tr was obtained by fitting of the exponential part of the 
distribution for each year. The values of Tr are shown 
in Fig. El The plots for the 1980s and 1990s are shifted 
vertically for clarity. We observe that the data points in 
the lower-income part of the distribution collapse on the 
same exponential curve for all years. This demonstrates 
that the shape of the income distribution for the lower 
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FIG. 8 Mam panel: Lorenz plots for income distribu- 
tion in 1983 and 2000. The data p oints are from the 
IRS (jStrudler. Petska. and Petskal . l2003l ). and the theoret- 
ical curves represent Eq. 1)20^ with the parameter / de- 
duced from Fig. [Tj Inset: Th e closed circles are the IRS 
data (,Str udler. Petska. and Petska . .2003 ) for the Gini coef- 
ficient G, and the op en circles show the theoret ical formula 
G = {1 + /)/2. From lSilva and Yakovenkcl (|2005l ). 



class is extremely stable and does not change in time, de- 
spite gradual increase in the average income in nominal 
dollars. This observation suggests that the lower-class 
distribution is in statistical "thermal" equilibrium. 

On the other hand, as Fig. [7] shows, income distribu- 
tion of the upper class does not rescal e and significantly 
changes in time. ISilva and Yakovenkol (|2005l ) found that 
the exponent a of the power law C(r) cx 1 /r" decreased 
from 1.8 in 1983 to 1.4 in 2000. This means that the up- 
per tail became "fatter" . Another useful parameter is the 
total income of the upper class as the fraction / of the to- 
tal income in the system. T he fraction / increased from 
4% in 1983 to 20% in 2000 (ISilva and Yakovenkol 120051) . 
However, in year 2001, a increased and / decreased, in- 
dicating that the upper tail was reduced after the stock 
market crash at that time. These results indicate that 
the upper tail is highly dynamical and not stationary. It 
tends to swell during the stock market boom and shrink 
during the bust. Similar results were found for Japan 
(lAovama et aLLl2003l:[Fuiiwara et aTl . l2003l:ISoumal . l200ll 
I2OO2D . 



Although relative income inequality within the lower 
class remains stable, the overall income inequality in the 
USA has increased significantly as a result of the tremen- 
dous growth of the income of the upper class. This is 
illustrated by the Lorenz curve an d the Gini coeffi cient 
shown in Fig. [5] The Lorenz curve ( KakwaiiH Il980[ ) is a 
standard way of representing income distribution in the 
economic literature. It is defined in terms of two coordi- 
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FIG. 9 Historical evolution of the parameters (r), Tr, and the 
fraction of income / going t o the upper tail, as defined in Eq. 
(pTjl . From lBaneriS (|2oS)- 



nates x(r) and y{r) depending on a parameter r: 

/;r'P(r')rf/ 



a;(r) 



P(r')dr', y(r) - 



(19) 



The horizontal coordinate x{r) is the fraction of the pop- 
ulation with income below r, and the vertical coordinate 
y(r) is the fraction of the income this population accounts 
for. As r changes from to oo, a; and y change from to 
1 and parametrically define a curve in the {x, y) plane. 
Fig. [5] shows the data points for the Lorenz 

curves in 19 83 and 2000, as com puted by 

the IRS (IStrudler. Petska. and Petskal . l2003f) . 
iDragulescu and Yakovenkon2001a[ ) analvticallv derived 
the Lorenz curve formula y = x + {1 — x) ln(l — x) for 
a purely exponential distribution P{r) = cexp(— r/T^). 
This formula is shown by the upper curve in Fig. [8] and 
describes the 1983 data reasonably well. However, for 
year 2000, it is essential to take into account the fraction 
/ of income in the upper tail, which modifies for the 



20051) 



Lorenz formula as follows I Dragulescu and Yakovenko . 
2003 : ISilva and Yakovenkol 20051 : lYakovenko and Silva . 



y^{l- f)[x+{l-x) ln(l - x)] + f Q{x - 1). (20) 



The last term in Eq. ((20)) represent the vertical jump of 
the Lorenz curve at a; = 1, where a small percentage of 
population in the upper class accounts for a substantial 
fraction / of the total income. The lower curve in Fig. [5] 
shows that Eq. (|20| fits the 2000 data very well. 

The deviation of the Lorenz curve from the straight 
diagonal line in Fig. [5] is a certain measure of income 
inequality. Indeed, if everybody had the same income, 
the Lorenz curve would be the diagonal line, because 
the fraction of income would be proportional to the frac- 
tion of the population. The standard measure of income 
inequality is the Gini coefficient < G < 1, which is 
defined as the area between the Lorenz curve and the 



diagonal line, divided by the area of th e triangle be- 
neath the diagonal line ( Kakwaiil Il980[ ). Time evo- 



lution of the Gini coefficient^, as co mputed by the IRS 
(Strudler. Petska. and Petska. 2003), is shown in the in- 



set of Fig. [H Dragulescu and Yakovenkol ( 2001al) derived 
analytically the result that G = 1/2 for a purely exponen- 
tial distribution. In the first approximation, the values 
of G shown in the inset of Fig. [5] are indeed close to the 
theoretical value 1/2. If we take into account the upper 
tail using Eq. (|^D|) , the formula for the Gini coefficie nt be- 
comes G = {l + f)/2 (ISilva and Yakovenkol [20051) . The 
inset in Fig. [5] shows that this formula gives a very good 
fit to the IRS data for the 1990s using the values of / 
deduced from Fig. H The values G < 1/2 in the 1980s 
cannot be captured by this formula, because the Lorenz 
data points are slightly above the theoretical curve for 
1983 in Fig. [S] Overall, we observe that income inequal- 
ity has been increasing for the last 20 years, because of 
swelling of the Pareto tail, but decreased in 2001 after 
the stock market crash. 

It is easy to show that the parameter / in Eq. (PO)) and 
in Fig. [5] is given by 



/ = 



(r) 



(r) 



(21) 



where (r) is the average income of the whole popula- 
tion, and the temperature is the average income in 
the exponential part of the distribution. Eq. (PT|) gives 
a well-defined measure of the deviation of the actual in- 
come distribution from the exponential one and, thus, 
of the fatness of the upper tail. Fig. [9] shows historical 
evolution of the parameters (r), Tr, and / given by Eq. 
((2T|) .^^ We observe that Tr has been increasing, essen- 
tially, monotonously (most of this increase is inflation). 
In contrast, (r) had sharp peaks in 2000 and 2006 coin- 
ciding with the speculative bubbles in financial markets. 
The fraction /, which characterizes income inequality, 
has been increasing for the last 20 years and reached 
maxima of 20% in the years 2000 and 2006 with a sharp 
drop in between. We conclude that the speculative bub- 
bles greatly increase the fraction of income going to the 
upper tail, but do not change income distribution of the 
lower class. When the bubbles inevitably collapse, in- 
come inequality reduces. 

Thus far we discussed the distribution of individual in- 
come. An interesting related question is the distribution 
P2 (r) of family income r = ri + r2, where ri and r2 are 
the incomes of spouses. If the individual incomes are 
distributed exponentially P{r) oc exp(— r/T^), then 



^2(0 = / dr'P{r')P{r -r') 
Jo 



crexp(-r/rr), (22) 



A similar plot was constructed by ISilva and Yakovenkol l|2005t ) 
for an earlier historical dataset. 
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United States, Bureau of Census data for 1996 




Annual family income, k$ 

FIG. 10 Histogram: Probability distribution of family in- 
come for families with two adults (US Census Bureau data). 
Solid li ne: Fit to Eq. (|22|l . From IPragulescu and Yakovenkol 
(|2001al 1. 



PSID data for families, 1 999 

100^ ^ ' ' 




Labor income of one earner, k$ 

FIG. 11 Scatter plot of the spouses' incomes (ri,r2) and 
(r2,ri) based on the da ta from the Panel Study of Income 
Dynamics (PSID). From [PragulesciT and Yakovenkol ||2003D . 



United States, Bureau of Census data for 1947-1994 



where c is a normalization constant. Fig. [TO] shows 
that Eq. is in good agreement with the family 

i ncome distribution data from the US Census Bureau 
(jPragulescu and Yakove'^ l2001al) . In Eq. dH]), we as- 
sumed that incomes of spouses are uncorrelatcd. This 
simple approximation is indeed supported by the scat- 
ter plot of incomes of spouses shown in Fig. [TT] Each 
family is represented in this plot by two points (ri,r2) 
and (r2 , ri ) for symmetry. We observe that the den- 
sity of points is approximately constant along the lines 
of constant family income ri + r2 = const, which indi- 
cates that incomes of spouses are approximately uncor- 
related. There is no significant clustering of points along 
the diagonal ri — r2, i.e., no strong positive correlation 
of spouses' incomes. 

The Gini coefficient for the family income dis- 
tribution (1^ was a.nalytica lly calculated by 
IPragules cu gmd Yakovenkol (|2001at ) as G = 3/8 = 37.5%. 



Fig. [T2] shows the Lorenz quintiles and the Gini co- 
efficient for 1947-1994 plotted from the US Census 
Bureau data ( Dragulescu and Yakovenkol l2001al ). The 
solid line, representing the Lorenz curve calculated from 
Eq. (im, is in good agreement with the data. The 
systematic deviation for the top 5% of earners results 
from the upper tail, which has a less pronounced effect 
on family income than on individual income, because 
of income averaging in the family. The Gini coefficient, 
shown in the inset of Fig. [121 is close to the calculated 
value of 37.5%. Notice that income distribution is 
very stable for a long per iod o f time, which was also 
recognized by economists ( Levvl . Il987t ). Moreover, the 
average G for the developed capitalist countries of North 
America and western Europe, as determined by the 
World Bank, is also close to the ca lculated value 37.5% 
( Dragulescu and Yakovenkol . |2003| ). However, within 




10 20 30 40 50 60 70 80 90 100% 
Cumulative percent of famiiies 

FIG. 12 Mam panel: Lorenz plot for family income calculated 
from Eq. (|22p . compared with the US Census data points. 
Inset: The US Census data points for the Gini coefficient 
for families, comp ared with the theoretically calculate d value 
3/8=37.5%. From lDraeulescu and Yakovenkol (|2001al ). 



this average, nations or groups of nations may have quite 
different Gini coefficients that persist over time due 
to specific historica l, pofi tical, or social circumstances 
(jRosser and Rosserj . |2004| ). The Nordic economies, with 
their famously redistributive welfare states, have G 
in the mid-20%, while many of the Latin American 
countries have G over 50%, reflecting entrenched social 
patterns inherited from the colonial era. 

Income distribution has been examined in 
econophys ics papers for different countries: 



Jagan (lAovama et all 
l2005t iFuiiwara etM 



Ferrero 



2004 . 



iNirei and Souma 



18 



2007 



2005), 



2005al : 



2007), 



ISoumal . I2OOIL 

Ger many 



l2002t iSouma and Nirei 
(IClementi and Gallegat 



the 



2005a 



2007 



IClementi. Gallegati. and Kaniadakis , 
UK (Clementi and Gallegati, 



Clementi. Gallep;ati. and Kaniadakii 



Ferrero 



2004 



200s 



Richmond. Hutzler . Coelho. and Repetowiczl._ 



2006 



2006 ). Italy^ (Clementi, Di Matt co. and Gallegati 



Clementi and Gallegati! 



Clementi. Gallegati. and Ka niadakid. 



2005hl: at^^o_ 



the USA 
iRawlings et al\ . 
Australia 



(Clementi and Gallegat il. 
2004 . India dsinhal. 




2006 



(iBaneriee. Yakovenko. and Di Matteo 



Clementi. Di Matteo. and Gallegatil . 12006 ; 
I)i Matteo. Aste. and Hvdd . 1200 4*). and New Zealand 
iFerrero . 2004. 2005ir The distributions are qualitatively 
similar to the results presented in this section. The 
upper tail follows a power law and comprises a small 
fraction of population. To fit the lower part of the 
distribution, different papers used the exponential. 
Gamma, and log-normal distributions. Unfortunately, 
income distribution is often reported by statistical 
agencies for households, so it is difficult to differentiate 
between one-earner and two-earner income distributions. 
Some papers used interpolating functions with different 
asymptotic behavior for low and high incomes, such as 
the Tsallis function (,Ferrero . 2005.) and the Ka n iadaki s 
function ( Clementi. Gallegati. and Kaniadakid . |2007| ) . 
However, the transition between the lower and upper 
classes is not smooth for the US data shown in Figs. 
[S] and [71 so such functions would not be useful in 
this case. The special case is income distribution in 
Argentina during the economic crisis, which shows a 
time- dependent bimodal shape with two peaks ( Ferrerol 
I2OO5D . 



B. Theoretical models of Income distribution 

Having examined the empirical data on income distri- 
bution, let us now discuss theoretical models. Income 
Ti is the influx of money per unit time to an agent i. 
If the money balance mi is analogous to energy, then 
the income Vi would be analogous to power, which is the 
energy flux per unit time. So, one should conceptually 
distinguish between the distributions of money and in- 
come. While money is regularly transferred from one 
agent to another in pairwise transactions, it is not typi- 
cal for agents to trade portions of their income. Never- 
theless, indirect transfer of income may occur when one 
employee is promoted and another demoted while the to- 
tal annual budget is fixed, or when one company gets a 
contract whereas another one loses it, etc. A reasonable 
approach, which has a long tradition in the economic 
literature ( Champernownd . Il953l : ICibratl . Il93ll : iKaleckil 
I1945D . is to treat individual income r as a stochastic pro- 
cess and study its probability distribution. In general, 
one can study a Markov process generated by a matrix of 



transitions from one income to another. In the case where 
the income r changes by a small amount Ar over a time 
period At, the Markov process can be treated as income 
diffusion. Then one can a pply the gener al Fokker-Planck 
equation ( Lifshitz and Pitaevskiii . [19811 ) to describe evo- 
lution in time t of the income distribution function P(r, t) 
(ISilva and YakovenkoL [2005h 



dP 
dt 



d_ 

dr 



AP 



d{BP) 



dr 



A = - 



(Ar) 
At 



B = 



2At 



(23) 

The coefficients A and B in Eq. ([23|l are determined by 
the first and second moments of income changes per unit 
time. The stationary solution dtP = of Eq. ((23)) obeys 
the following equation with the general solution 



d{BP) 
dr 



-AP, P{r) 



B{r) 



exp 



A{r') 
B{r') 



dr' 



(24) 

For the lower part of the distribution, it is reasonable 
to assume that Ar is independent of r, i.e., the changes 
in income are independent of inc ome itself. This pro- 
cess i s called the additive diffusion ( Silva and Yakovenkol 
I2OO5I ). In this case, the coefficients in Eq. ([23l) are the 
constants and B^. Then Eq. ([M)) gives the expo- 
nential distribution P{r) cx exp(— r/T^) with the effec- 
tive income temperature T,. B^/Af)}'^ The coincidence 
of this result with the Boltzmann-Gibbs exponential law 
([1]) and ([7]) is not accidental. Indeed, instead of con- 
sidering pairwise interaction between particles, one can 
derive Eq. ([T]) by considering energy transfers between 
a particle and a big reservoir, as long as the transfer 
process is "additive" and d oes not involve a Maxwel l- 
demon-like discrimination ('Basu a nd Mohantvl . l2008f) . 
Although money and income are different concepts, they 
may have similar distributions, because they are gov- 
erned by simila r ma thematical principles . It w as shown 
explicitly by [ Cord icr. Parcschi, and Toscani ([2 005); 



iDragulescu and Yakovenkol (|2000( ): ISlanhial (|2004D that 
the models of pairwise money transfer can be described 
in a certain limit by the Fokker-Planck equation. 

On the other hand, for the upper tail of income dis- 
tribution, it is reasonable to expect that Ar oc r, i.e., 
income changes are proportional to incom e itself. This i s 
known as the proportionality principle of lGibra"3 ( 193ll ). 
and the process is called the multiplicative diffusion 
(ISilva and Yakovenkol I200I . \ n this case, A = ar and 
B — br^, and Eq. ([24]) gives the power-law distribution 
P{r) oc l/r°+i with a = l + a/b. 

Generally, the lower-class income comes from wages 
and salaries, where the additive process is appropriate, 
whereas the upper-class income comes from bonuses, 
investments, and capital gains, calculated i n percent- 
ages, where the multiplicative process applies ( Milakov"i3, 



Notice that a meaningful stationary solution l|24^ requires that 
A>0, i.e., (Ar) < 0. 
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[2005') . However, the additive and multiplicative pro- 
cesses may coexist. An employee may receive a cost-of- 
living raise calculated in percentages (the multiplicative 
process) and a merit raise calculated in dollars (the addi- 
tive process) . Assuming that these processes are uncorre- 
lated, we have A = Ao+ar and B = Bo+br^ = 5(rg+r^), 
where Tq = Bo /b. Substituting these expressions into Eq. 



we find 



P(r) 



g-(ro/Tr) arctan(r/ro) 

[l-h(r/ro)2]i+°/26 



(25) 



The distribution (j25p interpolates between the exponen- 
tial law for low r and the power law for high r, be- 
cause either the additive or the multiplicative process 
dominates in the corresponding limit. The crossover be- 
tween the two regimes takes place at r = tq, where 
the additive and multiplicative contributions to B are 
equal. The distribution (P5)) has three parameters: the 
"income temperature" Tr = Aq/Bq, the Pareto expo- 
nent a = 1 + a/b, and the crossover income tq. It is 
a minimal model that captures the salient features of 
the empirical i n come distribution. Eq. (^5]) was obtained 
by lYakovenkol (|2009f) . and a more general formula for 
correlate d additive and mult i plicati ve processes was de- 
rived by iFiaschi and MarsilH ( 20091 ) for a sophisticated 
economic model. Fits of the IRS dat a using Eg. (1^ 
are shown in Fig. [12] reproduced from iBanerie^ (|2008f ). 
A mathematically similar, bu t more economica l ly or i- 



ented, model was propo sed by iNirei and" Soumal ((2003); 
ISouma and Nbdl (I2005D . where labor income and assets 
accumulation are described by the additive and multi- 
plicative processes correspondingly. A general stochastic 
process with additi ve and multiplic a tive n oise was stud- 



ied numerically by iTakavasu et Ml I|199<D . but the sta- 



tionary distribution was not derived analytically. A sim- 
ilar pro c ess w ith discrete time increments was studied by 
iKestenI (I1973D . Besides economic applications, Eq. ([^5)) 
may be also useful for general stochastic processes with 
additive and multiplicative components. 

To verify the multiplicative and additive hypotheses 
empirically, it is necessary to have data on income mo- 
bility, i.e., the income changes Ar of the same people 
from one year to another. The distribution of income 
changes P{Ar\r) conditional on income r is generally not 
available publicly, although it can be reconstructed by 
researchers at the tax agencies. Nevertheless, the mul- 
tiplicative hypoth esis for the up p er cla s s was quantita- 
tively verified by lAovama et al. I (|2003t) : iFuiiwara et al\ 
(|2003h for Japan, where such data for the top taxpayers 
are publicly available. 

The power-law distribution is meaningful only when 
it is limited to high enough incomes r > tq. If all 
incomes r from to oo follow a purely multiplicative 
process [Aq = and Bq = 0), then one can change 
to a logarithmic variable x — ln(r'/r:,,) in Eq. (|23p and 
show that it gives a Gaussian time-dependent distribu- 
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Annual Income/T 



FIG. 13 Fits of the IRS data for income distribution using 
Eq. (|25p . Pl ots for di f ferent years are shifted vertically for 
clarity. From iBaneried l|2008l '). 



tion Pt{x) (X exp( 
distribution for r, 



-X /2a t) for x, i.e., the log-normal 
also known as the Gibrat distribu- 



tion (i Gibratl . [ 193ll ). However, the width of this distri- 
bution increases in time, so the distribution is not sta- 
tionary. This was pointed out by iKaleckil (jl945[) a long 
time ago, but the log-normal distribution is still widely 
used for fitting income distribution, despite this funda- 
men tal logical flaw in its ju stification. In the classic pa- 
per, [Chamgernown^ (1953) showed that a multiplicative 
process gives a stationary power-law distribution when 
a boundary condition is imposed at tq ^ 0. Later, this 
result was rediscovered by econophvsicists (lLevv . p003t 
iLevv and Solomonl . Il996l : ISornette and Contl . Il997l ). In 
Eq. ((25|) . the exponential distribution of the lower class 
effectively provides such a boundary condition for the 
power law of the upper class. Notice also that Eq. ([^5]) 
reduces to Eq. (|18p in the limit tq — > with Bq = 0, but 
Ao^O. 

There are alternative approaches to income distribu- 
tion in econo mic literature. One of them, proposed by 
iLvdalll ( 1959t ). involves social hierarchy. Groups of peo- 
ple have leaders, which have leaders of the higher order, 
and so on. The number of people decreases geometri- 
cally (exponentially) with the increase in the hierarchical 
level. If individual income increases by a certain factor 
(i.e., multiplicatively) when moving to the next hierar- 
chical level, th en income distribution follows a power law 
(ILvdalll . [1951 . However, this original argument of Ly- 
dall can be easily modified to produce the exponential 
distribution. If individual income increases by a certain 
amount, i.e., income increases linearly with the hierarchi- 



20 



cal level, then income distribution is exponential. The 
latter process seems to be more realistic for moderate 
annual incomes below 100 k$ . A similar scenario is the 
Bernoulli trials ("Felleil. f 19661 ). where individuals have a 
constant probability of increasing their income by a fixed 
amount. We see that the deterministic hierarchical mod- 
els and the stochastic models of additive and multiplica- 
tive income mobility represent essentially the same ideas. 



V. CONCLUSIONS 

The "invasion" of physicists into economics and fi- 
nance at the turn of the millennium is a fascinating 
phenomenon. It generated a lively public debate about 
the role and future perspectives of econophysics, covering 
both theoretical and empirical issues. The econophysi- 
cist Joseph McCauley proclaimed that "Econophysics 
will displace economics in both the universities and 
boardrooms, simply because what is taught in economics 
classes doesn't work" ( Balll . 20061) . Althoug h there is 
some truth in his arguments (lMcCaulevl . [200l . one may 



consider a less radical scenario. Econophysics may be- 
come a branch of economics, in the same way as game 
theory, psychological economics, and now agent-based 
modeling became branches of economics. These branches 
have their own interests, methods, philosophy, and jour- 
nals. When infusion of new ideas from a different field 
happens, the main contribution often consists not in 
answering old questions, but in raising new questions. 
Much of the misunderstanding between economists and 
physicists happens not because they are getting different 
answers, but because they are answering different ques- 
tions. 

The subject of income and wealth distributions and 
social inequality was very popular at the turn of another 
century and is associated with the names of Pareto, 
Lorenz, Gini, Gibrat, and Champcrnowne, among others. 
Following the work by Pareto, attention of researchers 
was primarily focused on the power laws. However, 
when physicists took a fresh look at the empirical data, 
they found a different, exponential law for the lower part 
of the distribution. Demonstration of the ubiquitous 
nature of the exponential distribution for money, wealth, 
and income is one of the new contributions produced by 
econophysics.^^ The motivation, of course, came from 
the Boltzmann-Gibbs distribution in physics. Further 
studies revealed a more detailed picture of the two-class 



See, for exa mple. lAnglirJ ll2005l'): iBalll ll200d): ICarbone et al\ 

t; [Gallegati. ^een~Lux. and OrmerodI l|2006t) : 

j200q ); lLux| _ |2005|, |2009l '): iMcCaulevI 12009): 
ond, Chakrabarti, Chattcricc. and Angle (200a); iRosserl 
(2006, 2008b ); Stauffcr (2004 ); Yakovcnka (2009) . 
The exponential distribution is also u biquitous in the probabil- 
ity distributions of financ i al returns |Kleinert and Chen| . |2007 ; 
McCauley and Gunaratnd . l2003l ; ISilva. Prange. and Yakovenkc , 



distribution in a society. Although social classes have 
been known in political economy since Karl Marx, real- 
ization that they are described by simple mathematical 
distributions is quite new. Very i nteresting work wa s 
done by the computer scientist Ian i Wright! ( 20051 [20091 ) . 
who demonstrated emergence of two classes in an agent- 
based simulation of initially equal agents. This work 
has been further developed in the upcoming book by 
ICottrell. Cockshott. Michaelson. Wright, and Yakovenkol 



20041) and the growth rates of firms. 



(|2009D . mtegrating economics, computer science, and 
physics. 

Econophysics may be also useful for teaching of sta- 
tistical physics. If nothing else, it helps to clarify the 
foundations of statistical physics by applying it to non- 
traditional objects. Practitioners of statistical physics 
know very well that the major fascinating attraction of 
this field is the enormous breadth of its applications. 
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